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I. INTRODUCTION 



Semiclassical methods provide a very useful tool for 
the study of nonperturbative processes. Tunneling phe- 
nomena represent one of the most notable cases where 
semiclassical techniques can be used to obtain other- 
wise unattainable information on transition probabilities. 
There are many cases where the application of semiclassi- 
cal methods reduces to having to solve the classical equa- 
tions of motion in the Euclidean time domain, namely 
with t analytically continued to it. A famous example is 
WKB approximation to tunneling in quantum mechan- 
ics with one degree of freedom. In systems with multiple 
degrees of freedom there are other, more complex situa- 
tions, however, where the evolution along the imaginary 
time axis cannot be treated separately from the real time 
evolution which precedes and follows it. This happens, 
for example, in the study of tunneling in high-energy col- 
lisions in field theory [lj, |2| , where one must consider a 
system with definite particle number in the initial state, 
or in the study of chemical reactions where the initial 
system is in a definite quantum state 0. Then, in the 
search for complex-time solutions of the classical equa- 
tions of motion that dominate tunneling, one encounters 
a novel phenomenon which, if not properly handled, may 
make it impossible to obtain the appropriate solutions. 
We ran into this kind of problems in the study of an inter- 
esting class of field theoretical processes, where tunneling 
occurs between different topological sectors and is accom- 
panied by baryon number violation pi 01 However, 
the specifics of the model are immaterial for the features 
of tunneling which we plan to describe in this paper. In- 
deed, the fact that, with a field theoretical system, one 
is dealing with an infinite number of degrees of freedom 
is also irrelevant for the appearance of the phenomenon 
that we alluded to above: this phenomenon emerges as 
soon as one deals with more than one degree of freedom 
and in quite a range of applications of semiclassical tech- 
niques. For this reason we decided to illustrate it in this 
paper in the context of a simple quantum mechanical 
system with two degrees of freedom. The problems one 
encounters in the study of the field theoretical system and 
their solution will be described in a separate paper 0. 
Our hope is that the discussion presented here may help 
to shed some light on how semiclassical methods can be 
effectively applied in situations which a priory may seem 
intractable. 

Our considerations are of quite general validity, but, 
for purposes of clarity, it is useful to present them in 
the context of a definite model. We will take this to be 
the model of Ref. 0, || , namely a system formed by two 
particles of identical mass m, moving in one dimension 
and bound by a harmonic oscillator potential of frequency 
u> (Fig.^). One of the particles interacts with a repulsive 
potential barrier. This potential barrier is assumed to be 
high and wide, while the spacing between the oscillator 
levels is much smaller than the barrier height Vq. Thus, 
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FIG. 1: Oscillator hitting a potential barrier, with which only 
the "dark" particle interacts. 

the Hamiltonian of the model is 

where the conditions on the oscillator frequency and po- 
tential barrier are 

hw < Vo, (2) 
a > h/y/mVo . 

The problem is to calculate the transition probability 
through the barrier for a given incoming bound state; 
since the variables do not separate, this is certainly a 
non-trivial problem. Quantum-mechanically the incom- 
ing state is fully characterized by its total energy E and 
oscillator excitation number N. For E < Vq the transi- 
tion through the barrier can happen only through tunnel- 
ing. However, the transition may be classically forbidden, 
and occur through tunneling, also if E > Vq. A detailed 
study of the model (Ref. y| ) shows that one can dis- 
tinguish two regions in the E — N plane, as illustrated 
in Fig. |21 region A where transition across the barrier is 
classically forbidden, and region B where it is classically 
allowed. (Since Hui(N — |) cannot exceed E, the further 
region with A > E/huj- \ ~ E/h ui is trivially excluded.) 
The vertex S of B marks the lowest energy at which one 
can have a classically allowed transition, characterized by 
E = Es = Vq and some definite initial excitation number 
N = Ns- Strictly speaking, in order to have a classically 
allowed transition the energy must be (infinitesimally) 
larger than Es. For E = Es there exists a static, unsta- 
ble solution to the classical equations of motion, where 
the system just sits on top of the barrier: x\ — x 2 = 0. 
In the study of field theoretical tunneling processes this 
solution has been called the sphaleron [10| and, for sake 
of terminology, we will use the term also in this article, 
but no special meaning should be attributed here to it: 
for us it will be just a name. The line Epi in Fig.|2|marks 
a special class of tunneling solutions to the equations of 
motion, about which more will be said later. 

In the region A one would like to calculate the 
exponent F(E, N) in the transition probability T w 
exp{—F/K). The boundary between regions A and B is 
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FIG. 2: Allowed and forbidden regions in the model with two 
degrees of freedom. 



given by the line E (N), where F(E, N) = 0. F can be 
calculated by the method introduced in Ref. which 
will be recapitulated later. (See also the earlier work 
by Miller Q for an elegant approach to the calculation 
of the transition probability, based on canonical trans- 
formations.) The method of Ref. leads to having 
to solve the classical equations of motion, after analytic 
continuation to complex time and complex phase space, 
subject to suitable boundary conditions. In the tunneling 
solutions, the system evolves along the real time direction 
first, then undergoes the under-barrier evolution along 
the imaginary time direction. For the solutions situated 
on the line Epj in Fig. |5J the coordinates of the particles 
are real throughout the evolution, but this is a very spe- 
cial case. For general tunneling solutions the coordinates 
and the corresponding momenta, have to be continued 
to complex values. This is a defining feature of pro- 
cesses with more than one degree of freedom: in systems 
where a single particle moves along a line and tunnels 
under a potential barrier, all tunneling solutions would 
be real and behave like the processes on the line Epj. In 
systems with multiple degrees of freedom tunneling solu- 
tions, i.e. evolutions which obey the classical (analytically 
continued) equations of motion and obey the appropri- 
ate boundary conditions, can be rather readily found for 
E < Es- As the energy approaches the barrier height, 
however, something special happens: in the tunneling so- 
lutions the system lingers for progressively longer times 
on top of the barrier and the boundary conditions become 
truly asymptotic in time. Rigorously, finite time evolu- 
tions that satisfy the boundary conditions cease to exist. 
Here lies the subtlety in the application of semiclassical 
methods to such processes and the root of the difficulty 
in finding solutions that can be used throughout the do- 
main of classically forbidden transitions (i.e. the whole of 
region A). As we plan to show in this paper, the difficulty 
can be overcome by employing a suitable regularization, 
and one can find physically relevant solutions which map 
the entire region of classically forbidden processes and 
also merge continuously with classically allowed evolu- 
tions as one approaches the boundary Eq(N). The new 
solutions we find all describe tunneling onto the top of the 
potential barrier, or onto long-living excitations above 



the top. Notice that the problem is special to transitions 
which have an energy larger than the barrier height (i.e. 
larger than the sphaleron energy) and thus would never 
be encountered in systems with a single degree of free- 
dom, the systems normally used in textbooks to illustrate 
the use of semiclassical techniques, where all processes 
with energy larger than the barrier height are classically 
allowed. Still, even in simple systems with a single degree 
of freedom, one can see a sign of the subtleties that one 
can, and does encounter in more complex systems. For 
this reason, and also in order to put our whole formu- 
lation on a clear footing, we will start by re-examining 
the semiclassical description of tunneling in the case of 
a single particle moving on a line. This will be done in 
Sec. El 

In Sec. IIIII we consider the quantum mechanical model 
with two degrees of freedom. In particular, in Sec. IIII Al 
we introduce the semiclassical technique for the calcula- 
tion of the tunneling exponent. Then we examine the 
classical over-barrier solutions and find the boundary of 
the classically allowed region Eq(N) in Sec. IIII Bl In 
Sec. IIII CI we present a straightforward application of the 
semiclassical technique, outlined in Sec. IIII Al and find 
that it ceases to produce the relevant classical solutions 
in a part of the tunneling region A, roughly speaking, 
at E > Es- In Sec. II VI we introduce our regularization 
technique and show that it indeed enables one to find 
the tunneling classical solutions in the entire region A 
(Sec. lIV A*j l. We check our method against the numerical 
solution of the full Schrodinger equation in Sec. II V fil In 
Sec. lIV O we show how our regularization technique may 
be used to join smoothly the classically allowed and clas- 
sically forbidden regions. Our conclusions are presented 
in Sec. E| 



II. QUANTUM MECHANICS OF ONE DEGREE 
OF FREEDOM 

A. Classical boundary value problem 

Let us review the derivation of the semiclassical expres- 
sion for the tunneling exponent in quantum mechanics of 
one degree of freedom. The result is well known in this 
case (WKB formula) and the derivation here may seem 
artificially complicated, but it can be generalized to the 
case of multiple degrees of freedom, while many of its 
general features can still be captured in the relatively 
simple case of one degree of freedom. 

Let us calculate the probability of tunneling through 
a potential barrier, located in a region near X = 0, from 
the asymptotic region X — > — oo. Let us use a system of 
units with h — 1 and set the mass of a particle equal to 1 
(both can be reintroduced in our formulas on dimensional 
grounds). The Hamiltonian we consider has the general 
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form 



H 



where the potential V(X) has the form 




(3) 



V(X) = ^U{VXX) 
X 

and rapidly decays as X — > ±00. Then the semiclassical 
limit of this model is obtained when A — > while the 
energy is kept parametrically large: E = E/X. 

The probability of tunneling from an incoming (X — > 
—00) state with energy E is 



+00 



T(E) = lim 

t f — ti — >oo 



dX", 



-iH(t f -ti) 



\E) 



(4) 



where it is implicit that the initial state \E) is a wave 
packet with very small spread in momentum P, centered 
at large negative X. The transition amplitude 

A fi = (X f \e- iA ^-^\X i ) 

and its complex conjugate A*,f have the standard path 
integral representation 



A 



ft 



dX] 
dX' 



x(t i )=x i 

X(t f ) = X f 



(5) 



-iS[X'] 



x'(t i )=x' i 
x'(t f )=x f 



where S[X] = J. f L(X,X)dt is the action of the model. 
Inserting factors J dXi \Xi)(Xi\ — 1 one re-expresses (@J 
in terms of these transition amplitudes and the initial 
state matrix elements 

B iv = (X t \E)(E\X'J 

in the following way, 



+00 

T(E)= lim / dX f ( dXidX'.AfiAlfBi; 

tf-U->ca J J 



(6) 



To obtain an integral representation of the initial-state 
matrix elements, let us rewrite Ba> as follows, 



B w = (Xi\P E \X<) 



(7) 



where Pe denotes the projector onto a state with total 
energy E. In the momentum basis this projector has an 
integral representation which follows from the integral 
representation of the <5-function: 



(Q\Pe\P) 



s/2E 



d£e 



-*(J3-^-)« 



S(Q - P) 



(8) 



where \Q), \P) are the momentum eigenstates. It is con- 
venient to introduce the notation 



T : 



Combining formulas (JSJ , © , , iJSJ, and evaluating the 
Gaussian integrals over P and Q, one obtains the path 
integral representation for the tunneling probability. Af- 
ter rescaling X — > X/yX, X' — > X' j\[X it takes the 
following form 



y dT J [dX dX' 



exp 



-~F[X,X';T] 



>, (9) 



where the integration over Xj, X| and Xf is implicit, and 
X': T] = -iS[X] + iS[X'\ 

-ET+ ( Xi ~ X tf (10) 

with JB = XE being the rescaled energy (independent of 
A in the semiclassical limit A — ► 0). In this paper we will 
omit the tilde over the rescaled parameters, where this 
does not cause confusion. 

At small A one applies the saddle point approxima- 
tion to the integral 10. We will not consider the pre- 
exponential factor in the probability, so we will be inter- 
ested in the saddle point value of F. The saddle point 
equations for X(t) and X'{t) are the classical equations 
of motion, 



ss(x) = o ss(x') = o 



sx 



5X> 



(11) 



while the variation with respect to the boundary values 
Xf, Xi, X[ and auxiliary variable T leads to the bound- 
ary conditions 



Pf 
P, 

E -- 



P' 

Pi = 



Xi — X; 



iT 
2T 2 



(12a) 
(12b) 

(12c) 



where P = X. Note also, that one has (see Eq. |J5J) 

X{t f ) = X'(t f ) = X f , (12d) 

and that Eq. (|12cl) may be written in the form E = Pf /2. 
The latter equation requires that the energy of the solu- 
tion is equal to E. One more observation is that 



X' = X* 



(13) 



because X' corresponds to the saddle point of the com- 
plex conjugate amplitude. Equations l|12a|) and (|12d|) 



FIG. 3: Contour in the complex time plane. 



simply mean that the solutions X(t) and X'(t) coincide 
as t — > +00. On the other hand, Eq. I|12b|l states that 
the solutions X(t) and X'(t) are different at initial time 
t — *■ —00, if the saddle point value of T is not zero. There 
is no contradiction: for energy E larger than the barrier 
height, a real classical trajectory running from X — > —00 
to X — > +00 exists; for this solution T — and X co- 
incides with X' at any moment of time. On the other 
hand, if E is smaller than the barrier height, then there 
is no such classical solution in real time, but the solution 
may exist on a contour in the complex time plane, which 
wraps around some branch point (as shown in Fig. 
The solution continued along this contour is in general 
complex at large negative times (region A' in Fig. 0). 
Then equations l|12b|) , i|12c|) and guarantee that the 
saddle point value of T is real. 

The boundary conditions may be significantly simpli- 
fied by shifting the initial time from the real axis. As 
Ret — * — 00, the potential is irrelevant, and any solution 
describes free motion with momentum Pj = \f2E, 

X(t) = X, + P(t - U) , 

x'(t) = [x(t)r . 

Let us choose the initial time in such a way that Imt = 
T/2, i.e. consider the boundary value problem Ulljl . 112|) 
on the contour ABCD of Fig. El Then Eq. i(T2bT) takes a 
simple form 



1 1 
t = t + — , t = real — > —00 , 

which simply means that the coordinate is real on the 
part AB of the time contour. Equation l|13|) implies that 
the solution is real in the asymptotic region D, i.e. at t — 
real — > +00. Hence, instead of Eq. (|12d() one has reality 
conditions in the asymptotics A and D. After obtaining 
the solution to the classical boundary value problem, one 
inserts it to the right hand side of Eq. (|10|l and obtains 

F{E) = 2ImS- ET , (14) 

where S is the classical action on the contour ABCD. 

In quantum mechanics of one degree of freedom, the 
contour ABCD may be chosen in such a way that the 
points B and C are the turning points of the solution. 



FIG. 4: Classical solution describing tunneling. 

Then the solution is real also at the part BC of the con- 
tour. Indeed, a real solution at the part BC of the con- 
tour oscillates in the upside-down potential, T/2 equals 
to half-period of oscillations, and the points B and C are 
the two different turning points, X = (see Fig.0J). The 
continuation of this solution, according to the equation of 
motion, from the point C to the positive real times corre- 
sponds to the real-time motion, from the point Xc with 
zero initial velocity, towards X — > +00; the coordinate 
X(t) stays real on the part CD of the contour. Likewise, 
the continuation back in time from the point C leads to 
real solution in the part AB of the contour. In this way 
the reality conditions are satisfied at A and D. The only 
contribution to F comes from the Euclidean part of the 
contour, and one can check that the expression (|14|) re- 
duces to 

F{E)=2 J y/2(V(X) - E) dX , (15) 

X B 

which is the standard WKB result. 

To summarize, in quantum mechanics of one degree 
of freedom, the semiclassical calculation of the tunnel- 
ing exponent may be performed by solving the classical 
equation of motion on the contour ABCD in complex 
time plane, with the conditions that the solution be real 
in the asymptotic past (region A), and asymptotic future 
(region D). The relevant solutions tend to X — ► —00 and 
X — > +00 in regions A and D, respectively. The auxiliary 
parameter T is related to the energy of incoming state by 
the requirement that the energy of the classical solution 
equals to E. The exponent for the transition probability 
is given by Eq. I|14f> . which, in fact, coincides with the 
WKB result JTSJ. 

B. Joining the allowed and forbidden regions 

The solutions appropriate in the classically forbidden 
and classically allowed regions apparently belong to quite 
different branches. It is clear from Fig. 0] that as the en- 
ergy approaches the height of the barrier Vq from below, 
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the amplitude of the oscillations in the upside-down po- 
tential decreases, while the period T tends to a finite 
value determined by the curvature of the potential at its 
maximum. On the other hand, the solutions for E > Vq 
(classically allowed region) always run along the real time 
axis, so the parameter T is always zero. Hence, the rele- 
vant solutions do not merge at E = Vo, and T(E) has a 
discontinuity at E — Vq. We find it instructive to develop 
a regularization technique that removes this discontinu- 
ity and allows for smooth transition through the point 
E = Vq. The reason is that in quantum mechanics of 
multiple degrees of freedom, similar points exist not only 
at the boundary of the classically allowed region, but also 
well inside the forbidden region (but still at E > E$, see 
Introduction and Sec. IIII B"|) . 

To illustrate the situation, let us consider an exactly 
solvable model with 

U(X) = \r— . 

cosh 2 X 

The solution of the classical boundary value problem 
for energy smaller than the barrier height, E < Vq = 1, 
is 



E 



1 -E 



sinh A = - cosh ( V2E(t - t ) 



(16) 



This solution has branch points at t J , where the potential 
becomes singular (cosh A = 0): 



if — to — 



1 



±log 



E 



VT^E 



7T 

l 2 



(17) 



The requirement that the solution is real on the part AB 
of the contour determines the imaginary part of to, 



Imto = T/2 



(18) 



To fix the translational invariance along real time, one 
imposes a constraint on Re to- The solutions that are 
real at imaginary (Euclidean) time are obtained for 

Re t = . 

Later on, another choice will be convenient, for which the 
real parts of the position of the left singularities (t ~ ) are 
held at a certain negative value — c. This gives 



Re t 



1 



log 



i + Ve 



(19) 



We require X to be real on the CD part of the contour (in 
general, only asymptotically), and in this way determine 
the value of T/2. In the classically forbidden case, E < 1, 
we are interested in two branches of solutions. The first 
branch has T/2 = 0, i.e. the contour ABCD reduces to 
real axis; the classical solutions on this branch approach 
the barrier from the left and bounce back to X — > — oo. 



This branch is not relevant for tunneling. The tunneling 
branch 1 has 



T 
~2 



'2E 



(20) 



Equations l|18l) . 119|) and (|20|) determine the tunneling 
solution completely, while the transmission exponent is 
obtained by plugging this solution into Eq. 1141) . 

It is straightforward to obtain the classical solutions at 
E>1, 



E-l 



sinh A = sinh y\ 



to) 



(21) 



The role of the two branches is now reversed. The branch 
with T/2 = describes the over-barrier transitions, i.e., 
it is the relevant one. The branch with T/2 = -k /\/2E 
still exists, but it has X — > — oo at both part A and 
part D of the contour, and thus is irrelevant. We note in 
passing that all these solutions are real on the CD part of 
the contour not only asymptotically, but at finite times 
also. 

It is seen that the transition from the tunneling to 
the over-barrier solutions is indeed not smooth: T has a 
discontinuity at the energy equal to barrier height E — 
Vq = 1. The energy E = Vo = 1 is special also in the 
following respect. As the energy tends to Vq from below, 
the relevant classical solutions (with T/2 = n/^/2E) stay 
longer and longer near the top of the barrier. In the 
limit E — > 1 the time the solution spends near X = 
is infinitely long. This is reflected also in Eq. (|1 T|) : the 
distance between the branch points 



2 i 1 

lOg ; 

E 



E 



E 



(22) 



tends to infinity as E — *■ 1. The same situation oc- 
curs when the energy approaches the height of the bar- 
rier from above: the classical over-barrier solutions get 
"stuck" near X = 0. In a certain sense, all these solu- 
tions tend, as E — > Vo = 1, to the unstable static solution 
X(t) = — the sphaleron — describing a particle that sits 
on the top of the potential barrier. The latter property 
suggests a method to deal with this discontinuous behav- 
ior of the solutions with respect to energy. The idea is to 
regularize the classical boundary value problem in such 
a way that the classical solution X(t) = be unreach- 
able. We implement this idea by formally changing the 
potential 



U{X) 



'-U(X) , 



(23) 



1 There is in fact an infinite set of tunneling solutions with T/2 = 
(2n + l)w/V2E and reflecting solutions with T/2 = 2nir/V2E, 
but the imaginary part of their action is larger than that of the 
solutions we discuss, so their contribution to tunneling is expo- 
nentially small. 
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which leads to a corresponding change of the classical 
equations of motion. Here e is a real regularization pa- 
rameter, the smallest parameter in the model. At the 
end of the calculations one takes the limit e — * 0. 

We do not change the boundary conditions in our clas- 
sical problem, i.e., we still require that X(t) be real in the 
asymptotic future on the real time axis and that X(t') 
be real as t' — > — oo on the part A of the contour ABCD. 
Then the conserved energy is real. The sphaleron solu- 
tion X(t) = has now complex energy (because the po- 
tential is complex). Hence, the solutions to our classical 
boundary value problem necessarily avoid the sphaleron, 
and one may expect that the solutions behave smoothly 
in energy. 

We will have to say more about our regularization in 
Sec. lIVI here we show that it indeed does the job of con- 
necting the over-barrier solutions smoothly to the tun- 
neling ones. 

The general solution to the regularized problem can be 
obtained from (|16[) . 



E 



E 



sinhX = -cosh [V2E(t-t )) , 



where to is the integration constant. The value of Imio is 
fixed by the requirement that Im A = at positive time 
t — > +oo, 



T 

imin = 



1 



=arg[e 



E}. 



The condition that the coordinate A is real on the initial 
part AB of the contour gives the relation between T and 
E, 



T 
2 



1 



{7r + arg (e _le - E)} 



(24) 



For e = and E < 1, the result T/2 = it/\/2E is repro- 
duced. 

Let us analyze what happens in the regularized case 
in the vicinity of the would-be special value of energy, 
E = 1. It is clear from Eq. JUJ) that T is now a smooth 
function of E. Away from E = 1, Eq. I|24|l can be written 
as follows, 



T 
~2 



*2E{E - 1) 



forbidden region, 1 — E 3> e 
allowed region, E — 1 ^> e. 



(25) 



Deep enough in the forbidden region, when l-£>£, the 
argument in equation (|24() is nearly zero and we return 
to the original result H20(l . When E crosses the region 
of size of order e around E = 1, the argument rapidly 
changes from O(e) to — n, so that T/2 changes from tt/v2 
to nearly zero. Thus, at E > 1 we arrive to a solution 
which is very close to the classical over-barrier transition, 
and the contour is also very close to the real axis. This 
is shown in Figs. 03 We conclude that at small but 





FIG. 5: The contour in complex time plane for e / 0, £ < 1 
(upper figure), E = 1 (the middle figure) and E > 1 (lower 
figure). Arrows show the motion of the contour and branch 
points as energy increases. 




FIG. 6: T as function of E for zero and nonzero e, thick and 
thin lines, respectively. 



finite e, the classically allowed and classically forbidden 
regions merge smoothly. 

At E < 1, the limit e — > is straightforward. For 
E > 1 a somewhat more careful analysis of the limit 
e — ► is needed. From Eq. I|25|l one observes that the 
limit e — > with constant finite T < tt\/2 leads to so- 
lutions with E = 1. Classical over-barrier solutions of 
the original problem with E > E$ = 1 are obtained 
in the limit e — > provided that T also tends to zero 
while t = T/e is kept finite. Different energies corre- 
spond to different values of r. And that is what one 
expects — classical over-barrier transitions are described 
by the solutions on the contour with T = 0. 

Let us now discuss the behavior of the branch points 
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in complex time plane at small but finite e and varying 
energy. A particularly interesting question is whether or 
not the singularities cross our time contour ABCD. The 
positions of the singularities in the regularized problem 
are (we set Reio = here, see discussion after Eq. iJTSJl) 



1 f 7T 

- < inn + i h i ar 



g (e~ i€ / 2 ± y/E 



±ln 



- E 



The behavior of the real parts of the "right" (t+) and 
"left" (t~) singularities is simple — as the energy ap- 
proaches the barrier height, the singularities move to- 
wards large times, reaching a maximum distance from the 
origin equal to ln(2/e)/(2\/2) at E = 1, and then move 
back. This again shows that our regularization makes 
the solutions continuous in energy. The behavior of the 
imaginary parts of the positions of the "right" singular- 
ities, Im(i+), is also simple: they scale with energy like 
n /2\/2E and stay essentially constant as energy crosses 
the value E = 1. The imaginary parts of the positions of 
the "left" singularities, on the other hand, rapidly move 
down near E = 1, together with the contour height T/2 
(Fig. 0. So, the singularities do not cross the contour. 

Thus, by introducing small parameter e, Eq. (|23[) . we 
were able to regularize the system in such a way that the 
transition from classically forbidden region to the allowed 
one is rapid but smooth, and in the limit e — > only the 
physically relevant branches of solutions are selected. 

Although we have already studied all physically inter- 
esting solutions in non-regularized problem, let us ana- 
lyze what happens in the limit e — > to the regularized 
solutions with < T < ttv2 and E = 1. At this point it 
is convenient to make the choice l|19f) . Then in the limit 
e — > one obtains a family of solutions with intermediate 
T, 



sinh X 



3 -v^(t-ii+c) 



(26) 



These solutions are no longer real at any finite time, but 
they are still real asymptotically, as t — > +oo, where they 
approach the sphaleron, X = 0. For all these solutions 
the suppression exponent I|1U|) is zero. 

Thus there are three different branches of solutions 
merging at each of the bifurcation points E = 1, T/2 = 
and E = l, T/2 = tt/\/2. This is shown in Fig.0 We will 
see that a somewhat similar situation occurs in the case 
of multiple degrees of freedom, but there the branch anal- 
ogous to the branch is not degenerate in energy and 
is precisely the branch of relevant solutions in a certain 
region of parameters. The generalization of our regular- 
ization technique will enable us to select the physically 
relevant branches of solutions. 



' — — A: tunneling 


> — — - — _C: reflection 


- 


E: tunneling 




sphaleron 


D: reflection 


B: over-barrier 







FIG. 7: The dependence of the period T on energy E. Five 
branches of solutions are shown. A: represents tunneling so- 
lutions, Eq. \2\)\\ . in the forbidden region (E < Vo = 1); B: 
classically allowed transitions, T — 0; (C and D): unphysi- 
cal solutions undergoing reflection, in allowed and forbidden 
regions respectively; E new solutions 



III. QUANTUM MECHANICS OF TWO 
DEGREES OF FREEDOM 



The situation we discuss in this and the following sec- 
tion is a transition through a potential barrier of an os- 
cillator whose spacing between the levels is much smaller 
than the barrier height. This model was already de- 
scribed in the Introduction, see Q for the Hamiltonian. 

As in the one dimensional case, we choose units with 
H = 1, m = 1. The oscillator frequency ui is of order 1. 
The system is semiclassical, i.e. conditions © are satis- 
fied, if one chooses a = 1/V2A, Vo = 1/A, where A is 
a small parameter. Note that A plays a role similar to 
H, determining the magnitude of quantum fluctuations. 
At the classical level, this parameter is irrelevant: after 
rescaling the variables 2 X\ — > X\j\[\, x-i — > xi/\J\, the 
small parameter enters only through the overall multi- 
plicative factor 1/A in the Hamiltonian. Therefore, the 
semiclassical technique can be developed as an asymp- 
totic expansion in A. 

The properties of the system described by the above 
Hamiltonian are made clearer by replacing the variables 
x\, X2 with the center-of-mass coordinate X = [x\ + 
x<i)/\pl and the relative oscillator coordinate y = (xi — 
X2) / V2. In terms of the latter variables, the Hamiltonian 
takes the form 



H = & + — + —V 2 + - e - x i x +v) 2 / 2 . 



(27) 



2 To keep notations simple, we use the same symbols xi, X2 for 
the rescaled variables. 
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The interaction potential 



V iM EE lo • 2 



vanishes in the asymptotic regions X — > ±00 and de- 
scribes a potential barrier between these regions. At 
X — > ±00 the Hamiltonian l|27|) corresponds to an os- 
cillator of frequency u) moving along the center-of-mass 
coordinate X. The oscillator asymptotic state may be 
characterized by its excitation number N and total en- 



ergy E 



P'x 



+ u>(N + 1/2). We are interested in the 



transmissions through the potential barrier of the oscil- 
lator with given initial values of E and N . 



A. T/0 boundary value problem 

The probability of tunneling from a state with fixed 
initial energy E and oscillator excitation number N from 
the asymptotic region X — > —00 to any state in the other 
asymptotic region X — » +00 takes the following form: 



T(E, N) 



lim > 

-ti->oo ^— ' 

/ 



(28) 



where it is implicit that the initial and final states have 
support only well outside the range of the potential, with 
X < and X > 0, respectively. Semiclassical meth- 
ods are applicable when the initial energy and excitation 
number are parametrically large, 

E = E/X : N = N/\, 

where E and N are held constant as A — > 0. The transi- 
tion probability has the exponential form 



T = Dc 



where D is a pre-exponential factor, which we do not 
consider in this paper. Our purpose is to calculate the 
leading semiclassical exponent F(E, N). The exponent 
for tunneling from the oscillator ground state [ll], fl2L [jjj 
may be obtained by taking F(E, N) in the limit N — > 0. 

The exponent F(E, N) is again related to solutions of 
the complexified classical boundary value problem. The 
derivation is similar to that given in Sec. Ill Al we present 
the details in App. El The outcome is as follows. In ad- 
dition to the Lagrange multiplier T one now has another 
Lagrange multiplier 9 which is related to the new pa- 
rameter iV characterizing the incoming state. We again 
rescale the variables, X — > X/y/X, y — > y/V~X, and omit 
tilde over the rescaled quantities E, N. As above, the 
boundary value problem is conveniently formulated on 
the contour ABCD in the complex time plane (see Fig-EJ , 
with the imaginary part of the initial time equal to T/2. 
The coordinates X(t), y(t) must satisfy the complexified 



equations of motion on the internal points of the contour, 
and be real in the asymptotic future (region D): 



5S 



5S 



SX(t) 8y(t) 

lmy(t) ->0, 
lmX(t) -> 0, 



= 0, 



as t — ► +00 



(29a) 



(29b) 



In the asymptotic past (region A of the contour, where 
t = t' + iT/2, t' is real negative) one can neglect the 
interaction potential t/; n t, and the oscillator decouples: 



1 



y 



(ue 



The boundary conditions in the asymptotic past, t' — > 
—00, are that the center-of-mass coordinate X must be 
real, while the complex amplitudes of the decoupled os- 
cillator must be linearly related, 



ImX -> 0, 

V — > e e u*, 



as t' 



(29c) 



The boundary conditions l|29b|) and (|29c|l make, in fact, 
eight real conditions (since, e.g., lm.X(t') — > implies 
that both Iml and Iml tend to zero), and completely 
determine a solution, up to time translation invariancc 
(see discussion in Adp. IA^I . 

It is shown in App. lAl that a solution to this boundary 
value problem is, in fact, an extremum of the functional 
(cf. dJ) 

F[X,y; X*,y*; T,6] = -iS{X,y] + iS[X* ,y*} 

- ET - N6 + Boundary terms . (30) 

The value of this functional at the extremum gives the 
exponent for the transition probability 

F(E, N) = 21mSo(T, 6) — ET — N6 , (31) 

where So is the action of the solution (integrated by parts, 
see App. EJ, 



Do 



dt 



1 d 2 X 1 d 2 y 1 2 2 



The values of the Lagrange multipliers T and 9 are re- 
lated to energy and excitation number as follows, 



E(T,{ 



d_ 

df 
d 



2Im5*o(T,l 



N(T,6) = — 2Im5 (T, 



(32) 
(33) 



Making use of Eq. Ij31|l , it is straightforward to check also 
the inverse Legendre transformation formulas, 



T{E, N) 
9(E, N) 



d_ 

"dE 
d 
~dN 



F(E, N) 
F(E, N) 



(34) 
(35) 
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One can also check that the right hand side of Eq. i|32|) 
coincides with the energy of the classical solution, while 
the right hand side of Eq. l|3"3"jl is equal to the classical 
counterpart of the occupation number, 

X 2 

E = — + ujN ; N = uv . (36) 

So, one may either search for the values of T and 9 that 
correspond to some given E and N, or, following a com- 
putationally simpler procedure, solve the boundary value 
problem (|29l) for given T and 9 and then find the corre- 
sponding values of E and N from Eq. H36(l . 

Let us discuss some subtle points of the boundary value 
problem l|29|) . First, one notices that the condition of 
asymptotic reality (|29b|) does not always coincide with 
the condition of reality at finite time. Of course, if the so- 
lution approaches the asymptotic region X — > +oo on the 
part CD of the contour, the asymptotic reality condition 
(|29bfl implies that the solution is real at any finite posi- 
tive t. Indeed, the oscillator decouples as X — > +oo, so 
the condition l|29bjl means that its phase and amplitude, 
as well as X(t), are real as t — > +oo. Due to equations of 
motion, X(t) and y(t) are real on the entire CD-part of 
the contour. This situation corresponds to the transition 
directly to the asymptotic region X — > +oo. However, 
the situation can be drastically different if the solution 
on the final part of the time contour remains in the inter- 
action region. For example, let us imagine that the solu- 
tion approaches the saddle point of the potential X = 0, 
y = as t — > +oo. Since one of the excitations about 
this point is unstable, there may exist solutions which 
approach this point exponentially along the unstable di- 
rection, i.e. X(t), y(t) oc e~ const t with possibly complex 
pre- factors. In this case the solution may be complex at 
any finite time, and become real only asymptotically, as 
t — > +oo. Such solution corresponds to tunneling to the 
saddle point of the barrier, after which the system rolls 
down classically towards X — > +oo (with probability of 
order 1, inessential for the tunneling exponent F). We 
will see in Sec. IIV Al that the situation of this sort in- 
deed takes place for some values of energy and excitation 
number. 

Second, since at large negative time (in the asymptotic 
region X — > — oo) the interaction potential disappears, 
it is straightforward to continue the asymptotics of the 
solution to the real time axis. For solutions satisfying 
(|29c|) this gives at large negative time 

y{t) = — L (ue-^e-^ 1 + u*e 9+ ^e iUJt ) , 
T 

lmX(t) = -i—px ■ 

We see that the dynamical coordinates on the negative 
side of the real time axis are generally complex. For so- 
lutions approaching the asymptotic region X — > +oo as 
t — > +oo (so that X and y are exactly real at finite t > 0), 
this means, in complete analogy to the one-dimensional 



case, that there should exist a branch point in the com- 
plex time plane: the contour in Fig. winds around this 
point and cannot be deformed to the real time axis. This 
argument does not work for solutions ending in the in- 
teraction region as t — > +oo, so branch points between 
the AB-part of the contour and the real time axis may 
be absent. We will see in Sec. IIV Al that this is indeed 
the case in our model for a certain range of E and N . 



B. Over-barrier transitions: boundary of the 
classically allowed region Eo(N) 

Before studying the exponentially suppressed transi- 
tions, let us consider the classically allowed ones. To this 
end, let us study the classical evolution, in which the 
system is initially located at large negative X and moves 
with positive center-of-mass velocity towards the asymp- 
totic region X — > +oo. The classical dynamics of the sys- 
tem is specified by four initial parameters. One of them 
(e.g., the initial center-of-mass coordinate) fixes the in- 
variance under time translation, while the other three are 
the total energy E, the initial excitation number of the 
y-oscillator, defined in classical theory as N = E osc /ui, 
and the initial oscillator phase ipi. 

The classically allowed region is the region of the E- 
N plane where over-barrier transitions are possible for 
some value(s) of the initial oscillator phase 3 . For given 
N, at large enough E the system can certainly evolve to 
the other side of the barrier. On the other hand, if E 
is smaller than the barrier height, the system definitely 
undergoes reflection. Thus, there exists some boundary 
energy Eq{N), such that classical transitions are possible 
for E > E (N), while for E < E (N) they do not occur 
for any initial phase tp.i. The line Eq(N) represents the 
boundary of the classically allowed region. We have cal- 
culated Eq(N) numerically for lu = 0.5: the result is 
shown in Fig. and reproduced in Fig. El 

An important point of the boundary Eq(N) corre- 
sponds to the static unstable classical solution X(t) = 
y(t) = 0, the sphaleron (cf. Refs. [H Q). It is 
the saddle point of the potential U(X, y) = cu 2 y 2 /2 + 
U- m t(X,y) and has exactly one unstable direction, the 
negative mode (see Fig. |SJ). The sphaleron energy E$ — 
C/(0,0) = 1 determines the minimum value of the func- 
tion Eq(N). Indeed, while classical over-barrier transi- 
tions with E < Eg are impossible, the over-barrier solu- 
tion with slightly higher energy can be obtained as fol- 
lows: at the point X = y = 0, one adds momentum along 
the negative mode, thus "pushing" the system towards 
X — > +oo. Continuing this solution backwards in time 
one finds that the system tends to X — > — oo for large 



3 Note that the corresponding classical solutions obey the bound- 
ary conditions I29bl . I29cl with T = 9 = 0, i.e., they are solu- 
tions to the boundary value problem 1291 . 
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negative time and has a certain oscillator excitation num- 
ber. Solutions with energy closer to the sphaleron energy 
correspond to smaller "push" and thus spend longer time 
near the sphaleron. In the limiting case when the energy 
is equal to Eg , the solution spends an infinite time in the 
vicinity of the sphaleron. This limiting case has a defi- 
nite initial excitation number Ng, so that Eq(Ns) = E$ 
(see Fig. [SJ). The value of iVg is unique because there 
is exactly one negative direction of the potential in the 
vicinity of the sphaleron. 

In complete analogy to the features of the over-barrier 
classical solutions near the sphaleron point (Eg, N$), one 
expects that as the values of E, N approach any other 
boundary point (Eq(N), TV), the corresponding over- 
barrier solutions will spend more and more time in the 
interaction region, where U ln t 7^ 0. This again follows 
from a continuity argument. Namely, let us first fix the 
initial and final times, U and tf. If in this time interval 
a solution with energy E\ evolves to the other side of 
the barrier and a solution with energy E 2 and the same 
oscillator excitation number is reflected back, there ex- 
ists an intermediate energy at which the solution ends 
up at t = tf in the interaction region. Taking the limit 
tf —>■ +00 and {E% — E2) —> 0, we obtain a point at the 
boundary of the classically allowed region and a solution 
tending asymptotically to some unstable time dependent 
solution that spends infinite time in the interaction re- 
gion. We call the latter solution excited sphaleron; it de- 
scribes some (in general nonlinear) oscillations above the 
sphaleron along the stable direction in coordinate space. 
Therefore, every point of the boundary (Eq (N) , N) cor- 
responds to some excited sphaleron. 

We display in Fig. [S] a solution, found numerically in 
our model, that tends to an excited sphaleron. We see 
that this trajectory is, roughly speaking, orthogonal to 
the unstable direction at the saddle point (X = 0, y = 0). 



C. Suppressed transitions: bifurcation line E\(N) 

Let us now turn to classically forbidden transitions, 
and consider the boundary value problem (|29|) . It is rel- 
atively straightforward to obtain numerically solutions 
for 9 = 0. The boundary conditions l|29bl) . (|29c[) in this 
case take the form of reality conditions in the asymptotic 
future and past. It can be shown (see Ref. ^5|) that the 
physically relevant solutions with 9 = are real on the 
entire contour ABCD of Fig.|21and describe nonlinear os- 
cillations in the upside-down potential on the Euclidean 
part BC of the contour. The period of the oscillations is 
equal to T, so that the points B and C are two different 
turning points, where X = y = 0. These real Euclidean 
solutions are called periodic instantons (see Ref. [l5|). A 
practical technique for obtaining these solutions numer- 
ically on the Euclidean part BC consists in minimizing 
the Euclidean action (for example with the method of 
conjugate gradients, see Ref. @ for details). The solu- 
tions on the entire contour are then obtained by solving 




y+X 



FIG. 8: The potential (dotted lines) in the vicinity of the 
sphaleron (X = 0, y = 0) (marked by the point), excited 
sphaleron (thick line) corresponding to the point (E, N) = 
(1.985,3.72) at the boundary of the classically allowed re- 
gion, and the trajectory of the solution which is close to this 
excited sphaleron (thin line). In this figure the asymptotic 
regions X — > ±00 are along the diagonal. 




0.5 E s = l 1.5 

E 



FIG. 9: The boundary of the classically allowed region Eo(N), 
the bifurcation line E\ (N) and the line of the periodic instan- 
tons Epi(N). 



numerically the Cauchy problem, forward in time along 
the line CD and backward in time along the line BA. 
Having the solution in asymptotic past (region A), one 
then calculates its energy and excitation number l|36|l . 
The solutions to this Cauchy problem are obviously real, 
so the boundary conditions (|29bp . (|29c|) are indeed sat- 
isfied with 9 = 0. It is worth noting that solutions with 
9 = arc similar to the ones in quantum mechanics of 
one degree of freedom, see Fig. 0J The line of periodic 
instantons in E-N plane in our model is shown in Fig. [5J 
again for uj = 0.5. 

Once the solutions with 9 = are found, it is nat- 
ural to try to cover the entire classically forbidden re- 
gion of the E — N plane with a deformation procedure, 
by moving in small steps in 9 (and T). The solution 
to the boundary value problem with (T + AT, 9 + A6>) 
may be obtained numerically, by applying an iteration 
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FIG. 10: The dependence of the tunneling coordinate X 
on time for two solutions with nearly the same energy and 
initial excitation number. The physical solution tunnels to 
the asymptotic region X — > +oo, while the unphysical one 
gets reflected back to X — > — oo. The physical solution 
has E = 1.028, N = 0.44, while the unphysical one has 
E = 1.034, N = 0.44. These two solutions are close to the 
point at the bifurcation line Ei(N = 0.44) = 1.031. 



technique, with the known solution at (T, 9) serving as 
the zeroth-order approximation 4 . Provided the solutions 
have correct "topology" at each step, i.e. X — > — oo and 
X — > +oo on parts A and D of the contour, the solutions 
obtained by this "walking" procedure are physically rel- 
evant. However, the "walking" method fails to produce 
relevant solution if there are bifurcation points, where 
the physical branch of solutions merges to an unphysical 
branch. As there are unphysical solutions close to physi- 
cal ones in the vicinity of a bifurcation point, one cannot 
"walk" near these points. 

In our model, the walking method produces correct 
solutions to the T/9 boundary value problem in a large 
region of the E — N plane, where E < Ei(N). How- 
ever, at sufficiently high energy E > E\(N), where 
Ei(N) > Es, the deformation procedure generates so- 
lutions that bounce back from the barrier (see Fig. I10[) . 
i.e. have the wrong "topology" . Clearly, the latter do not 
describe the tunneling transitions of interest. Therefore, 
if the semiclassical method is applicable at all in the re- 
gion Ei(N) < E < Eq(N), there exists another, physical 
branch of solutions. In that case the line Ei(N) is the 
bifurcation line, where the physical solutions "meet" the 
ones with wrong "topology" . Walking in small steps in 
9 and T is useless in the vicinity of this bifurcation line, 
and one needs to introduce some trick to find the relevant 
solutions beyond that line. The bifurcation line E\(N) 
for our quantum mechanical problem of two degrees of 
freedom with uj = 0.5 is shown in Fig. EI 

The loss of topology beyond a certain bifurcation line 
in E — N plane is by no means a property of our model 
only. This phenomenon has been observed in field theory 



4 In practice, the Newton-Raphson method is particularly conve- 
nient (see Refs. 0, Q, 0). 



models, in the context of both induced false vacuum de- 
cay and baryon-number violating transitions in gauge 
theory Q (in field-theoretic models, the parameter N is 
the number of incoming particles). In all cases, the loss 
of topology prevented one to compute the semiclassical 
exponent for the transition probability in an interesting 
region of relatively high energies, where the suppression 
of tunneling may be fairly weak. 

Coming back to the quantum mechanics of two de- 
grees of freedom, we point out that the properties of tun- 
neling solutions approaching the bifurcation line Ei(N) 
from the left, are in some sense similar to the proper- 
ties of tunneling solutions in one-dimensional quantum 
mechanics near the boundary of the classically allowed 
region. Again by continuity, these solutions of our two- 
dimensional model spend a long time in the interaction 
region; this time tends to infinity on the line Ei(N). 
Hence, on any point of this line, there is a solution that 
starts in the asymptotic region left of the barrier, and 
ends up on an excited sphaleron. We already mentioned 
in Sec. IIIlAl that such a behavior is indeed possible be- 
cause of the existence of an unstable direction near the 
(excited) sphaleron, even for complex initial data. We 
have suggested in Sec. Ill Bl a trick to deal with this 
situation — this is our regularization technique. 



IV. REGULARIZATION TECHNIQUE 

In this Section we further develop our regularization 
technique, and find the physically relevant solutions be- 
tween the lines Ei(N) and E (N). We will see that all 
solutions from the new branch (and not only on the lines 
Eq(N) and Ei(N)) correspond to tunneling onto the ex- 
cited sphaleron. These solutions would be very difficult, 
if at all possible, to obtain directly, by solving numeri- 
cally the non-regularized classical boundary value prob- 
lem (|29|) : they are complex at finite times, and become 
real only asymptotically as t — ► +oo, whereas numerical 
methods require working with finite time intervals. 

As an additional advantage, our regularization tech- 
nique enables one to obtain a family of the over-barrier 
solutions, that covers all the classically allowed region, in- 
cluding its boundary. This may be of interest in models 
with large number of degrees of freedom and in field the- 
ory, where finding the boundary of the classically allowed 
region by direct methods is difficult (see e.g., Ref. 16] for 
discussion of this point). 



A. Regularized problem in the classically forbidden 

region 

The main idea of our method is to "regularize" the 
equations of motion by adding a term proportional to a 
small parameter e so that configurations staying for an 
infinite time near the sphaleron no longer exist among 
the solutions of the T/9 boundary value problem. After 
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performing the regularization we explore all the classi- 
cally forbidden region without crossing the bifurcation 
line. Taking then the limit e-s Owe reconstruct the cor- 
rect values of F, E and N. 

When formulating the regularization technique it 
is more convenient to work with the functional 
F[X, y; X*, y*; T, 9] itself rather than with the equations 
of motion. We slightly modify the form of F, Eq. (|30|l . 
so that F is no longer extremized by configurations ap- 
proaching the excited sphalerons asymptotically. To 
achieve this, we add to the original functional l|30|) a 
new term of the form 2eT mt , where T ult estimates the 
time the solution "spends" in the interaction region. 
The parameter of regularization e is the smallest one 
in the problem, so any regular extremum of the func- 
tional F (the solution that spends finite time in the region 
C^int 7^ 0) changes only slightly after the regularization. 
At the same time, the excited sphaleron configuration 
has Tint = oo which leads infinite value of the regularized 
functional F e = F + 2eT lnt - Hence, the excited sphalerons 
are not stationary points of the regularized functional. 

For the problem at hand, £/i n t ~ 1 in the interaction 
region, and Tj nt can be defined as follows, 



dtU int {X,y)+ / dtU int (X*,y*) 



(37) 



We notice that Ti nt is real, and that the regularization is 
equivalent to the multiplication of the interaction poten- 
tial by a complex factor 



U int -* (1 - ie)U iDt = c- l£ U int + 0{e 2 ) 



(38) 



This results in the corresponding change of the clas- 
sical equations of motion, while the boundary condi- 
tions l|29bjl . (|29c|) remain unaltered. 

We still have to understand whether solutions with 
e ^ exist at all. The reason for the existence of such 
solutions is as follows. Let us consider a well-defined (for 
e > 0) matrix element 



lim 

-t 4 — > 



(-iH-eUi 



■>*)(*/-**) | .E, N) 



f 



where \E, N) denotes, as before, the incoming state with 
given energy and number of particles. The quantity T e 
has a well defined limit as e — > 0, equal to the tunneling 
probability (|28p. As the saddle point of the regularized 
functional F e = F + 2eT- mt gives the semiclassical expo- 
nent for the quantity T e , one expects that such saddle 
point indeed exists. 

Therefore, the regularized T/9 boundary value prob- 
lem is expected to have solutions necessarily spending 
finite time in the interaction region. By continuity, these 
solutions do not experience reflection from the barrier, if 
one makes use of the "walking" procedure starting from 
solutions with correct "topology" . The line Ex (N) is no 
longer a bifurcation line of the regularized system, so the 



walking procedure enables one to cover the entire forbid- 
den region. The semiclassical suppression factor of the 
original problem is recovered in the limit e — > 0. 

It is worth noting that the interaction time is Legendre 
conjugate to e, 



ld_ 



F e (E,N,e) 



(39) 



Useful equations can also be obtained from Eqs. (|32[1 and 



E = E\ t=0 + 2e^ + O{e) 
N = N\ c=a + 2e^ + 0(e) 



(40) 
(41) 



These equations show how the values of E and N change 
with e — away from the would-be bifurcation line E\{N) 
the values of E and N depend slightly on e, while at 
E ~ E\(N) where T- lnt takes large values (for small e), 
the dependence is strong and results in a considerable 
shift of the point (E,N). Eqs. gDJl, I0TJ may be 

used as a check of numerical calculations. 

We implemented the regularization procedure numer- 
ically. To solve the boundary value problem, we make 
use of the computational methods described in Ref. ^| . 
In our calculations we used w = 0.5. To obtain the 
semiclassical tunneling exponent in the region between 
the bifurcation line E±(N) and boundary of the allowed 
region, Eq(N), we began with a solution to the non- 
regularized problem deep in the forbidden region (i.e., 
at E < Ei(N)). For that value of T and 9 we increased 
the value of e from zero to a certain small positive num- 
ber. Then we changed T and 9 in small steps, keeping 
e finite, and found solutions to the regularized problem 
in the region E 1 (N) < E < E {N). These solutions 
had correct "topology" , i.e. they indeed ended up in the 
asymptotic region X — > +oo. Finally, we lowered e and 
extrapolated F, E and N to the limit e — > 0. 

The lines T = const for e = 0.002 are shown in Fig. II II 
We see that at e ^ these lines are smooth and extend 
from 9 = (the line of periodic instantons) to 9 = oo 
(the line N — 0). They indeed cover all the tunneling 
region without any irregularities at the bifurcation line 
Ei(N). 

The same lines T = const, but now in the limit e —> 
are shown in Fig. the lines 9 — const, again in the 
limit e — > 0, are presented in Fig. ^] It is seen that 
the functions T(E, N) and 9(E, N) are continuous, while 
their first derivatives have discontinuities at the bifur- 
cation line E±(N). As T and 9 are the derivatives of 
the tunneling exponent F(E, N), see Eqs. (j2J, (l3*5|l . the 
suppression exponent F itself has discontinuities of the 
second derivatives only. The lines F = const are shown 
in Fig.^] We see that the function F and its first deriva- 
tives are indeed continuous. 

Let us consider more carefully the solutions in the re- 
gion Ex(N) < E < Eq(N) which we obtain in the limit 
e — > 0. They belong to a new branch, and thus may 
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FIG. 11: Lines T = const for e = 0.002. Numbers near each 
line show the values of T. The boundary of the classically 
allowed region, the bifurcation line and the line of periodic 
instantons are shown by dashed lines. 
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FIG. 13: The lines 6 = const and d = const in the limit e — > 0. 
Numbers near the lines are the values of 9 and #. 




FIG. 14: The lines F — const at energies above Es = 1. 
Dashed line shows Ei(N). 



FIG. 12: The lines T = const (in the forbidden region) and 
r = const (in the allowed region), see Sec. II V CI in the limit 
e — > 0. Numbers near the lines are the values of T and r. 



exhibit new physical properties. Indeed, we found that, 
as the value of e decreases to zero, the solution at any 
point (E,N) with E 1 (N) < E < E (N) spends more 
and more time in the interaction region. The limiting 
solution corresponding to e = has infinite interaction 
time: in other words, it tends, as t — > +oo, to one of 
the excited sphalerons. The resulting physical picture is 
that at large enough energy (i.e., at E > Ei(N)), the 
system prefers to tunnel exactly onto an unstable classi- 
cal solution, excited sphaleron, that oscillates about the 
top of the potential barrier! To demonstrate this, we 
have plotted in Fig. ^1 the solution x{t) = (X(t),y(t)) 
at large times, after taking numerically the limit e —>■ 0. 
To understand this figure, one recalls that the potential 
near the sphaleron point X = y = has one positive 



mode and one negative mode. Namely, by introducing 
new coordinates c+, c_ , 



X 

V = - 
ctg2a 



cos a c + -\ 
- sin a c + 

' 2 ' 



sin a c_ 
h cos a C- 



one writes, in the vicinity of the sphaleron, 



where 



H = 1 + f 



uj± = — 1 



pi 



2 



-ci- 



^c 2 



Since the solutions to the T/8 boundary value problem 
are complex, the coordinates c + and c_ are complex too. 
We show in Fig. 1151 real and imaginary parts of c + and 
c_ at large real time t (part CD of the contour). Wc 
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FIG. 15: The large-time behavior of a solution with e = 
at (E = 1.05, N = 0.43). The coordinates X and y are 
decomposed in a basis of the eigenmodes near the sphaleron. 
Note that Imc + = 0. 



see that while Rec+ oscillates, the unstable coordinate 
c~ approaches asymptotically the sphaleron value: c_ — > 
as t — > +oo. The imaginary part of c_ is non-zero 
at any finite time. This is the reason for the failure of 
straightforward numerical methods in the region E > 
Ei(N): the solutions from the physical branch do not 
satisfy the conditions of reality at any large but finite 
final time. We have pointed out in Sec. lIII Al that this can 
happen only if the solution ends up near the sphaleron, 
which has a negative mode. 

It is interesting to compare the branches of solutions in 
the two- and one-dimensional models (the latter was dis- 
cussed in Sec.|nJ|. The dependence T(E) in these models 
is shown in Fig. ^| (for the two dimensional problem the 
graph corresponds to constant oscillator excitation num- 
ber N = 0.1). We see that the structure of the branches is 
similar, though in the one-dimensional case the solutions 
corresponding to "tunneling" onto the sphaleron are de- 
generate in energy and are not really useful for finding 
the transition exponent. In the two-dimensional model, 
on the other hand, the similar branch is the one relevant 
at Ex(N) < E < E (N). 




1.0 E,(N=0J) 



3 E (N=0.1) 1.5 1.6 



(a) 



"tu n tiding" 
spbuFnm 



(b) 



FIG. 16: Dependence of the parameter T — —dF/dE on 
energy for (a) the two dimensional model with fixed N = 0.1 
and (b)the one dimensional model. The lower picture is the 
same as Fig.0 Different lines correspond to different branches 
of classical solutions to T/6 boundary value problem. The 
branches labelled "reflection" have wrong "topology" . 



B. Regularization technique versus exact 
quantum— mechanical solution 

The quantum mechanics of two degrees of freedom is a 
convenient testing ground for checking the semiclassical 
methods |8|,|9j and, in particular, our regularization tech- 
nique. The solutions to the full stationary Schrodinger 
equation may be found numerically in this case at finite 
values of the semiclassical parameter A, and the results 
may be extrapolated to A — > 0. The suppression fac- 
tor may be then compared to the semiclassical result. 
We performed this check in the region E > E$ = 1, 
which is most interesting for our purposes. We applied 
the numerical techniques for solving the full stationary 
Schrodinger equation, which were developed in Ref. 8]; 
our results agree with Ref. y|] where comparison is possi- 
ble. The results of the full quantum mechanical calcula- 



tion of the suppression exponent F in the limit A->0 are 
represented by points in Fig. 1171 The lines in that figure 
represent the values, for constant N, of the semiclassical 
exponent F(E,N), which we obtained by making use of 
the regularization procedure and extrapolation to e — ► 0. 
We see that in the entire forbidden region (including the 
region E > Ei(N)) the semiclassical result for F coin- 
cides with the exact one. 



C. Classically allowed region 

Finally, let us show that our regularization procedure 
enables one to obtain a subset of classical over-barrier so- 
lutions, that exist at high enough energies. This subset 
is interesting, as it extends all the way to the boundary 
of the classically allowed region, E = Eq(N). In princi- 
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FIG. 17: The tunneling exponent F(E, N) in the region 
E > Es = 1. The lines show the semiclassical results 
while points represent the exact ones, obtained by solving 
the Schrodinger equation. The lines across the plot are the 
boundary of the classically allowed region Eo(N) and the bi- 
furcation line Ei(N). 



pie, finding this boundary is purely a problem of classical 
mechanics, and, indeed, in mechanics of two degrees of 
freedom one obtains this boundary numerically by solv- 
ing numerically the Cauchy problem for given E and iV 
and all possible oscillator phases, see Sec. IIII Bl How- 
ever, if the number of degrees of freedom is much larger, 
this classical problem becomes quite complicated, as one 
has to span a high-dimensional space of Cauchy data. 
As an example, a stochastic Monte Carlo technique was 
developed in Ref. 0] to deal with this problem in field 
theory context. The approach below may be viewed as 
an alternative to the Cauchy methods. 

First, let us recall that all classical over-barrier solu- 
tions with given energy and excitation number satisfy 
the T/9 boundary value problem with T = 0, 9 = 0. We 
cannot reach the allowed region of E — N plane without 
regularization, because we have to cross the line Eq(N) 
corresponding to the excited sphaleron configurations in 
the final state. However, the excited sphalerons no longer 
exist among the solutions of the regularized boundary 
value problem at any finite value of e. This suggests that 
the regularization enables one to enter the classically al- 
lowed region and, after taking an appropriate limit, ob- 
tain classical solutions with finite values of E, N. 

By definition, the classically allowed transitions have 
F = 0. Thus, one expects that in the allowed region, the 
regularized problem has the property that F e (E,N) = 
ef(E,N) + 0(e 2 ). In view of the inverse Legendre for- 
mulas <|34[1 . i|35|) the values of T and 9 must be of order e: 
T = er(E, N),6 = ed{E, N), where the quantities r and 
$ are related to the initial energy and excitation number 



(see Eq. l|52fl. l|55fll in the following way, 

E = lim —2 Im S e = ^2T int (r, i?) , (42) 
d d 

N = m 2lm Se = dd 2Tint(r ' ^ (43) 

and we have used the form (|38|l of the regularized po- 
tential. Therefore, one expects that one can invade the 
classically allowed region by taking a fairly sophisticated 
limit e — > with r = T/e = const, -d = 9/e = const. In 
the allowed region the parameters r and i? are analogous 
to T and 9. 

By solving the regularized T/9 boundary value prob- 
lem one constructs a single solution for given E and N. 
On the other hand, for e — there are more classical 
over-barrier solutions: they form a continuous family la- 
beled by the initial oscillator phase. Thus, after taking 
the limit e — > one obtains a subset of over-barrier solu- 
tions, which should therefore obey some additional con- 
straint. It is almost obvious, that this constraint is that 
the interaction time Tint, Eq. (|37|) . is minimal. This is 
shown in App. [5] 

The subset of classical over-barrier solutions obtained 
in the limit e — > of the regularized T/9 procedure ex- 
tends all the way to the boundary of the classically al- 
lowed region. Let us see what happens when one ap- 
proaches this boundary from the classically allowed side. 
At the boundary Eq(N), the unregularized solutions tend 
to excited sphalerons, so the interaction time Tj nt is in- 
finite. This is consistent with 142|) . (|43[l only if r and i9 
become infinite at the boundary. Thus, to obtain a point 
of the boundary one takes the further limit, 

(E (N),N)= lim (E(t,#),N(t,#)) . 

r j {}— const 
r — > + oo 

Different values of r/i? correspond to different points of 
the line Eq(N). In this way one finds the boundary of the 
classically allowed region without an initial-state simula- 
tion. 

We have checked this procedure numerically. The line 
t = const, $ = const, with varying e, is shown in Fig. 1181 
The limit e — » exists indeed — the values of E and N 
tend to the point (c) of the classically allowed region. 
The phase of the tunneling coordinate X{t) in complex 
time plane is shown in Fig. 1191 for the three points (a), 
(b) and (c) of the curve Fig. The branch points of 
the solution 5 , the cuts and the contour are clearly seen 
on these graphs. 

It is worth noting that the left branch points t*_ move 
down as T and 9 approach zero. Solutions to the right of 
the point (b) (along the curve of Fig, IT8*|) have left branch 



5 The phase of the tunneling coordinate turns by 7r around the 
branch point. The points where the phase of the tunneling coor- 
dinate turns by 2ir correspond to the zeroes of X(t). 
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E 

FIG. 18: The line r = T/e = const = 380, ■d = 9/e = const = 
130 connecting the points (a): T — 3.8, 9 = 1.3, e = 0.01 
and (c): T = 0, 9 = 0, e = 0. Point (b) in the middle has 
T = 1.82, 9 = 0.62, e = 0.0048. As the values of r and t? 
are quite large, the limiting point is close to the boundary of 
the classically allowed region. The reflection boundary and 
the boundary of the classically allowed region are shown by 
dashed lines. 

point in the lower complex half-plane, Imtl < 0. There- 
fore, the corresponding contour may be continuously de- 
formed to the real time axis. These solutions still satisfy 
the reality conditions asymptotically (see Fig. I15|l . but 
show nontrivial complex behavior at any finite time. 

Making use of the regularized T/9 procedure, one is 
able to approach the boundary of the classically allowed 
region from both sides. The points at this boundary may 
be obtained by taking the limits T — > 0, T/9 — const of 
the tunneling solutions and r — > +oo, t/& = const of 
the classically allowed ones. As t* = t/i? = T/9 by 
construction, the lines r* = const are continuous at the 
boundary Eq(N), though may have discontinuity of the 
derivatives. These lines are shown in Fig. |2JJ| We see 
that the variable r* can be used to parametrize the curve 
E (N). 

For completeness, the lines r = const and i? = const in 
the classically allowed region are presented in Figs. FLU ITU 



V. CONCLUSIONS AND DISCUSSION 

Though the semiclassical processes studied in this pa- 
per were dubbed "tunneling" , our results suggest some- 
what different interpretation at energies exceeding the 
minimum height of the barrier, Eg- At these energies, 
it is energetically allowed that the system jumps above 
the barrier, and restarts its classical evolution from the 
region near its saddle point. This is precisely what we 
have observed. From the physical viewpoint, this is not 
quite what is normally meant by "tunneling through a 
barrier" . Yet the transitions remain exponentially sup- 




(c) 



FIG. 19: The phase of the tunneling coordinate in complex 
time plane at the three points of the curve r = const, i? = 
const (see Fig. 1181 . The asymptotics X — > -co and X — * 
+oo correspond to arg(X) = 7r and correspondingly. The 
contour in the time plane is plotted with white line. 
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FIG. 20: Lines T/9 = const are continuously connected with 
the lines r/t? = const. The numbers near the lines give the 
values of T/9 = r/0. 
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pressed (until the energy reaches the boundary of the 
allowed region), but, intuitively, the reason is different: 
to jump above the barrier, the system has to undergo 
considerable rearrangement, unless the incoming state 
is chosen properly (i.e. unless the oscillator excitation 
number in our model is large enough). This rearrange- 
ment costs exponentially small probability factor. We 
note that similar exponential factor was argued to ap- 
pear in various field th eory processes with multi-particle 
final states [H E [3 H| . 

The T/9 boundary value problem, equipped with our 
regularization technique, enables one to deal with this 
situation, which, in the language of complex classical so- 
lutions, corresponds to a new physically relevant branch. 
The regularization procedure chooses automatically the 
correct branch, and thus provides an efficient way to cross 
the bifurcation lines where different branches merge. The 
phenomenon discussed in this paper appears to be quite 
general, as it exists in diverse quantum systems, from 
two-dimensional quantum mechanics to quantum field 
theory. The appropriate generalization of our regular- 
ized T/9 procedure is straightforward; we are going to 
report on this application to baryon-number violating 
processes in gauge theory in future publication 0. 
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APPENDIX A: T/9 BOUNDARY VALUE 
PROBLEM 



The semiclassical method for calculating the probabil- 
ity of tunneling from a state with a few parameters fixed 
was developed in (llE^Jril IT^ in context of field theoret- 
ical models and in 8, 9] in quantum mechanics. Here we 
outline the method, adapted to our model of two degrees 
of freedom. 



1. Path integral representation of the transition 
probability 

We begin with the path integral representation for 
the probability of tunneling from the asymptotic region 



X — > — oo through a potential barrier. Let the incoming 
state \E, N)s have fixed energy and oscillator excitation 
number, and has support only for X <gC 0, well outside 
the range of the potential barrier. The inclusive tunnel- 
ing probability for states of this type is given by 



{-^OO TOO 
/ dX f J dy f 
-oo 



(X f ,y f \e-^-^\E,N)\ J. (Al) 

This probability can be reexpressed in terms of the tran- 
sition amplitudes 

A fi = {X f , y f \e- lH ^-^\X u Vl ) (A2) 
and initial-state matrix elements 

B iV = (X u Vi \E, N)(E, N\Xl, y[) (A3) 
in the following way, 

{+oo 
J dX f J dXi dX[ 
-oo 

+ 00 v 

J dy l dy' l dy 1 A fl A*, f B lV \ . (A4) 

— oo 

The transition amplitude and its complex conjugate have 
the familiar path integral representation: 



•Afi = j [dx] 
A*, f = f [dx' 



iS[i 



(A5) 



iS[x'\ 



S'(t / ) = x / 



where x — (X, y) , and S is the action of the model. To 
obtain a similar representation for the initial-state matrix 
elements let us rewrite Bui as follows, 



B iv = (X t , yi \P E P N \X' u y[) , 



(A6) 



where Pm and Pg denote the projectors onto states with 
oscillator excitation number N and total energy E re- 
spectively. It is convenient to use the coherent state for- 
malism for the y-oscillator and choose the momentum 
basis for the X-coordinate. In this representation, the 
kernel of the projector operator PePn takes the form 



(q, b\P E P N \p, a) = 



1 



(2tt)2 
i 



d£ drj 



exp ( -iEt, - iNrj + -p 2 £ + e luJ ^ +l71 ba ) S(q - p) , 
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where \p, a) is the eigenstate of the center-of-mass mo- 
mentum px and y-oscillator annihilation operator a with 
eigenvalues p and a respectively. It is straightforward to 
express this matrix element in the coordinate represen- 
tation using the formulas 



(X\p) 



7T 



JpX 



'2it 



Evaluating the Gaussian integrals over a, b, p, q, we ob- 
tain 



/ 



Bii' — d£ dr] exp — iE£ - iNr] 



2 € 



A2 



Y Q — 2iu^—2irj 



(A7) 



where we omit the pre-exponential factor depending on 
77, £. For the subsequent formulation of the boundary 
value problem it is convenient to introduce the notation 



T = -ii 



= —17] 



Then, combining together the integral representations 
(|A7|) and (|A5|) and rescaling coordinates, energy and ex- 
citation number x — > x '/VA, E -> E/X, N -> N/X, we 
finally obtain 



{-t-zoo 
J dTd9 J [dxdx] 



exp{-jF[x, x'- T,ff\}}, (A8) 



where 



F[x, T, 6>] = y] + iS[X ', y'} 

-ET-N9 + B l (x l ,x' l ;T,e). (A9) 

Here the non-trivial initial term is 

"1 



2T 



2 _ e 2uiT+29 



(A10) 



In (|A8(I x and x' are independent integration variables, 
while x'j = Xf, see Eqs. (|A4ll . 



2. The boundary value problem 

For small A, the path integral (|A8|) is dominated by a 
stationary point of the functional F. Thus, to calculate 
the tunneling exponent, we extremize this functional with 
respect to all variables of integration: X(t), y(t), X'(t), 
y'(t), T, 9. Note that because of the limit tf — ti — * +00, 
the variation with respect to the initial and final values 
of coordinates leads to boundary conditions imposed at 
asymptotic t — > ±00, rather than at finite times U,tf. 
Note also that the stationary points may be complex. 

The variation of the functional l|A9(l with respect to 
the coordinates at intermediate times gives second order 
equations of motion, in general complexified, 



SS 



6S 



SS' 



SS' 



SX(t) Sy(t) 6X'(t) 5y'(t) 







(Alia) 



-00 are 



The boundary conditions at the final time t f 
obtained by extremization of F with respect to Xf = X'p 
yf = y'j. These are 



X 



f 



X 'f 



Vf 



y'f 



(Allb) 



It is convenient to write the conditions at the initial 
time (obtained by varying Xi, yi, X[, y[) in terms of 
the asymptotic quantities. At the initial moment of time 
U — > —00, the system moves in the region X — > — 00, well 
outside of the range of the potential barrier. Equations 
(|Alla() in this region describe free motion of decoupled 
oscillator, and the general solution takes the following 
form, 



X(t)=Xi+pi(t-ti), 



y(t) 



-iw(i-ti) 



while the solution for X'(t), y'(t) has similar form. For 
the moment, a and a are independent variables. The 
initial boundary conditions in terms of the asymptotic 
variables Xi, pi, a, a take the form: 



Pi = Pi 

a' + a' 
a + a = 



Xi - X[ 
iT ' 
^ T+e + a 
-uT-6 



(Allc) 



a'e" T+e 



The variation with respect to the Lagrange multipliers T 
and 9 gives the relation between the values of E, N and 
initial asymptotic variables (here we use the boundary 
conditions IjAllcIl ). 



E 



EL 
2 

N = ad 



uiN , 



(Alld) 



Equations (jAllafl - I jAlldjl constitute the complete set 
of saddle-point equations for the functional F. 
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The variables X 1 and y' originate from the conju- 
gate amplitude A*,* (see Eq. l|A5jll. which suggests that 
they are the complex conjugate to X, y. Indeed, the 
Ansatz X'(t) = X*(t), y'(t) = y*(t) is compatible with 
the boundary value problem IjAllfl . Then the Lagrange 
multipliers T, 9 are real, and the problem I|A11J| may 
be conveniently formulated at the contour ABCD in the 
complex time plane (see Fig- El ■ 

Now we have only two independent complex variables 
X(t) and y(t), that have to satisfy the classical equations 
of motion in the interior of the contour, 



SS 



ss 



6X(t) Sy(t) 



= 



(Al2a) 



The final boundary conditions (see Eq. i)Allb|0 become 
the conditions of the reality of the variables X(t) and 
y(t) at the asymptotic part D of the contour: 



ImXf — 0, Imy/ = 
ImXf = Imyy = , 



-oo . (A12b) 



The seemingly complicated initial conditions ijAllcj! sim- 
plify when written in terms of the time coordinate t' = 
t+iT/2 running along the part AB of the contour. Let us 
again write the asymptotics of a solution, but now along 
the initial part AB of the contour: 



X = X + Po (t' - U) , 

y = 



ue 



iu(t'-ti) , ve i<^(t' -U) 



In terms of Xq, yo, u and v, the boundary conditions 
(|Allc|) are 



Im Jq = 0, Impo = , 



(A12c) 



* 9 

u e 



Finally, we write Eqs. ijAlldfl in terms of the asymptotic 
variables along the initial part of the contour: 



E -- 
N 



Pi 
2 

: LUUV 



■ujN 



(A13) 



These equations determine the Lagrange multipliers T, 9 
in terms of E, N. Alternatively, we can solve the prob- 
lem (|A12|) for given values of T, 9 and find the values of 
E,N from Eqs. (|A13|) . what is more convenient compu- 
tationally. 

Given a solution to the problem (|A12(1 . the exponent 
F is the value of the functional l|A9(l at this saddle point. 
Explicitly 

F = -ET - N9 + 2ImS (T,9) , (A14) 
where So is the action of the system, integrated by parts 

1 d 2 



Sq 



dt 



1 d 2 

2d£ X ~2 y dT 2 



-y 2 -U int (X,y) 



Note that we did not make use of the constraints (|A13|I to 
obtain the formula l|A14[) . so we still can extremize (1 A 14(1 
with respect to T and 9 (see discussion in Sec. IIII A"j) . 

The classical problem (|A12|) is conveniently dubbed 
T/9 boundary value problem. Eqs. (|A12b|) and l|A12c|) 
imply eight real boundary conditions for two complex 
second-order differential equations lA12a|) . However, one 
of these real conditions is redundant: Eq. i|A12b|) im- 
plies that the (conserved) energy is real, so the condition 
Impo — > is automatically satisfied (note that the os- 
cillator energy E osc = louv = ue s uu* is real). On the 
other hand, the system (|A12(I is invariant under time 
translations along the real axis. This invariance may 
be fixed, e.g., by demanding that Re(X) takes a pre- 
scribed value at a prescribed large negative time t' (note 
that other ways may be used instead. In particular, for 
E < Ei(N) it is convenient to impose the constraint 
KeX(t = 0) = 0). Together with the latter requirement, 
we have exactly eight real boundary conditions for the 
system of two complexified (i.e. four real) second-order 
equations. 



APPENDIX B: A PROPERTY OF SOLUTIONS 
TO T/6 PROBLEM IN THE ALLOWED REGION 

For given E, N there is only one over-barrier classi- 
cal solution which is obtained in the limit e — * of the 
regularized T/9 procedure. To see what singles out this 
solution, let us analyze the regularized functional 



F e [q] =F[q}+2eT int [q], 



(Bl) 



where q denotes the variables x(t), x'(t) and T, 9 to- 
gether. The unregularized functional F has a valley of 
extrema q e (if) corresponding to different values of the ini- 
tial oscillator phase cp. Clearly, at small e the extremum 
of F t is close to a point in this valley with the phase 
extremizing T int [q e (if)] . 



dip 



T int [q e ^)}=0 



(B2) 



Hence, the solution qf of the regularized T/9 boundary 
value problem tends to the over-barrier classical solution, 
with Tint extremized with respect to the initial oscillator 
phase. 

Because U- ln t(x) > 0, T mt is a positive quantity with 
at least one minimum. In normal situation there is only 
one saddle point of F e , so by solving the T/9 boundary 
value problem one obtains the classical solution with the 
time of interaction minimized. 

The dependence of the interaction time on the initial 
oscillator phase for given values of energy and excitation 
number is shown in Fig. [2] We see that Tj nt indeed 
has one minimum, and the corresponding solution indeed 
coincides with the limit of the solution of the regularized 
T/9 boundary value problem (black point on the graph). 
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FIG. 21: The dependence of the interaction time Tint on the 
initial oscillator phase at fixed values of E and N, correspond- 
ing to the point (c) in Fig. 1191 The black point marks the 
initial oscillator phase and interaction time calculated on the 
solution (c). 
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